Thermodynamics of a Trapped Unitary Fermi Gas 
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We present the first model-independent comparison of recent measurements of the entropy and 
of the critical temperature of a unitary Fermi gas, performed by Luo et al., with the most complete 
results currently available from finite temperature Monte Carlo calculations. The measurement of 
the critical temperature in a cold fermionic atomic cloud is consistent with a value T c — 0.23(2)ej? 
in the bulk, as predicted by the present authors in their Monte Carlo calculations. 

PACS numbers: 03.75.Ss 
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The study of the properties of a Fermi gas in the 
unitary regime (when the s-wave scattering length a 
is large compared to the average interparticle separa- 
tion) emerged as one of the most fascinating theoreti- 
cal many-body problems since it was first formulated by 
G.F. Bertsch as the Many-Body X (MBX) challenge in 
1999 The experimental investigation of the uni- 

tary Fermi gas (UFG) began with its realization in cold 
atomic traps by O'Hara et al. at Duke University three 
years later Q. At unitarity (often referred to as "at res- 
onance"), when a — > 00, the properties of such a sys- 
tem are governed by deceptively simple laws. In par- 
ticular, the ground state energy per particle is given by 
E/N = 3£f£/5; where Ef — fr 2 k F /2m is the Fermi en- 
ergy of a noninteracting Fermi gas with the same number 
density n = N/V = k F /3ir 2 . The determination of the 
dimensionless constant £ was the subject of the MBX 
challenge and the best current accepted value was deter- 
mined a bit later through restricted/fixed node Monte 
Carlo (MC) calculations as £ = 0.42(1) 0,0, HQ- This 
value was confirmed by the zero temperature extrapola- 
tion of unrestricted MC calculations of Ref. [§], where 
£ = 0.44(3) was obtained. Theoretically, it was also 
found that this system is superfluid at low temperatures 
and the value of the pairing gap was estimated at zero 
temperature to be A = 0.504(24)e F 0, H, @]. (For lack 
of space we quote and comment here only on theoretical 
results obtained in controlled MC calculations, where the 
errors are typically, though not always, only of statistical 
origin. In all other theoretical approaches that we are 
aware of, the errors are essentially impossible to quan- 
tify due to the lack of any identifiable small parameter.) 
A number of finite temperature thermodynamic prop- 
erties of the homogeneous phase was determined as well 
d, 0, 03, 11 1 , even though there is still some disagreement 
concerning the exact value of the critical temperature T c , 
on which we shall comment later. The temperature de- 
pendence of the pairing gap has not been determined yet. 
On the experimental side there is a quite wide spread in 
values of the dimensionless parameter £ [11] determined 
in various experiments. However, the latest experiments 



seem to converge, possibly guided by the existence of 
firm theoretical results, to the expected value: 0.74(7) 
1, 0.51(4) [H, 0.32+2^ 0, 0.36(15) 0, 0.46(5) 0, 
0.45(5) [!3|, 0.41(15) 0. The measurements of the pair- 
ing gap are still in their infancy. Although it has been 
conclusively demonstrated that a UFG is superfluid at 
sufficiently low temperatures [l!| , the value of the pair- 
ing gap has only been determined so far in one experi- 
ment [20]. Moreover, the extracted value is significantly 
smaller than the theoretical value 0, d, 0| ■ Such a small 
pairing gap is inconsistent with the value of T c measured 
independently T3,21|. In Ref. [2~J] . T r - was found by first 



determining the entropy and the energy of such a system, 
a procedure which allows to establish an absolute temper- 
ature scale. In this experiment a number of properties of 
the atomic cloud were determined in the unitary regime. 
The cloud was then adiabatically brought to the BCS side 
of the Feshbach resonance, where kpa = —0.75 (here kp 
is the Fermi momentum corresponding to the central den- 
sity of the cloud). In this regime (ftp|a| < 1) one can use 
many-body perturbation theory [22j, together with the 
Local Density Approximation (LDA), to evaluate vari- 
ous cloud properties (see Ref. 0] for a comparison of 
MC results with perturbative many-body results). Mea- 
surements of the energy of the cloud can thus be related, 
using theory, to temperature and entropy. Since the en- 
tropy is conserved, by measuring the energy of the cloud 
one can determine the energy-entropy dependence in the 
unitary regime and its absolute temperature as well (from 
T = dE(S,N)/dS). Most of the atomic trapping po- 
tentials used in these experiments can be approximated 
rather well with harmonic potential wells. Such poten- 
tials can be shown to satisfy the virial theorem at uni- 
tarity, namely E{T,N) = 2N(U) = 3muj 2 (z 2 ) and 
therefore simply measuring the spatial shape of the cloud 
allows for a unique determination of the UFG energy at 
any temperature. One of the main goals of the present 
work is to present a critical analysis of the results of this 
experiment 2l|, in the light of available finite tempera- 
ture MC calculations. This work thus represents the first 
model-independent comparison of a full theory directly 



2 



with experiment. 

At unitarity (1/kpa — 0) the pressure of a homoge- 
neous UFG is determined by a universal function h T (z): 
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where T and /i are the temperature and the chemical po- 
tential, respectively. Remembering that the grand canon- 
ical potential is fi(V, T, fi) — —VV(T,fi) one can easily 
show that the energy of the system reads: E = 3W/2, 
where V is the volume of the system. Following a rea- 
soning 24] similar to that in Ref. [25| one can show that 
thermodynamic stability implies positivity hr(z) > 
and convexity h'^(z) > 0. In the high-temperature 
limit fi — » — oo and V(T,fi) tends from above to the 
free Fermi gas pressure. In the low-temperature limit 
V(T, fi) tends from above to P(Q, £,e F ) = 4/3e 5 J 2 £,/5. Sim- 
ilarly, at all temperatures the pressure calculated in the 
BCS/meanfield approximation will give a variational es- 
timate from below of V(T,fi). In Fig. Q] we illustrate 
these statements and plot the results of three finite tem- 
perature MC calculations 0, [l(| ■ It is not possible to 



extract the data for this plot from Ref. [llj. The results 
of Ref. [HI stand apart from the rest of the theoreti- 
cal and experimental results. While the results [|| agree 
with these bounds, three out of the six calculated points 
in Ref. (one of the lowest T/highest z and two at the 
highest T/lowest z's) slightly violate them. The point at 
fi/T « 3.24 is approximately where the authors of Ref. 
9] claim that the normal-superfluid phase transition oc- 
curs. 

We use our MC results [H (with a spatial lattice size 
8 x 8 x 8) to generate smooth interpolation formulas for 
the energy, chemical potential and entropy (see inset of 
Fig. [T]). Standard manipulations show that all the UFG 
thermodynamic potentials can be expressed in terms of a 
single function of one variable, a property known as uni- 
versality d,H,[2^|. This property was incorporated in our 
interpolation. At high temperatures we notice that our 
results smoothly approach the corresponding free Fermi 
gas results with some offsets for the energy, chemical po- 
tential and entropy [§]. 

At this point we assume that the LDA can be used 
to describe the properties of an atomic cloud in a trap. 
There has been no systematic study of the accuracy of 
LDA in the unitary regime. We can, however, easily es- 
timate the role of the gradient corrections for a noninter- 
acting Fermi gas in an anisotropic harmonic trap. Using 
methods described in 27J one can show that the ground 
state energy of a two-component fermion system in an 
anisotropic harmonic trap is given by: 
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where fl — {u x uj y uj z ) 1 l' i and w = {uj 2 . + uj 2 + u} 2 )/3Q. 



FIG. 1: MC data from Ref.[| (blue circles), Ref.0 (six black 
points) and Ref. [hj (black crosses). The four straight lines 
starting at the origin are the T — > limits of hr(z — ► oo) = 
2 2/5 z/£ 3/5 , where £ = 0.25(2) 0, f = 0.42 [1,0, £ = 0.59 for 
meanfield/BCS approximation and £ = 1 for the free Fermi 
gas model respectively. The two solid lines (red/lower and 
green/higher) correspond to hr(z) calculated in the free Fermi 
gas and the BCS/meanfield approximation hr(z) respectively. 
In the inset we show the fits to MC data [8|. 



In this formula the first term is the leading LDA con- 
tribution. Naively, one would expect the next to lead- 
ing order LDA correction to be proportional to N. It is 
a peculiarity of harmonic potentials, however, that the 
leading gradient corrections start at next order instead, 
namely at 0(N 2 ^ 3 ) [23] ■ At the same order one finds the 
so-called shell correction to the energy, given in this for- 
mula by HO{3N) 2 / 3 e sc {N)/8. The function e sc (N) has 
a vanishing average over particle number, is minimum 
when a shell is filled and maximum in the middle of a 
shell. This term depends strongly on the asymmetry 
of the harmonic potential and has maximum amplitude 
for spherical potentials. In such case the amplitude of 
e sc (N) is about 0.5, while it is significantly smaller for 
asymmetric wells. Finite temperatures [2tJ and pairing, 
even at unitarity [28j], have a smoothing effect on the 
shell energy correction term, but do not affect in any 
other major way the remaining leading gradient correc- 
tion term. At temperatures close to T c one would ex- 
pect the gradient corrections to play a noticeable role 
in the description of the pairing properties [2^]. How- 
ever, close to T c the pairing energy represents a relatively 
small contribution to the total energy, and so the errors 
in the total energy from using naive LDA around T c are 
likely to be small too. In a trap, the fraction of particles 
that are close to loosing superfluidity (namely for which 
x(r = T/ef{y)) < x c = 0.23(2),where x c is the criti- 
cal temperature in natural units Q) is also small. All in 
all, it appears that for the mostly-harmonic traps used in 
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typical experiments the role of the gradient corrections is 
relatively small and LDA is a reasonable approximation. 

In this approach, the grand canonical thermodynamic 
potential for a UFG confined by an external potential 
C/(r) is a functional of the local density n(r) given by 



n= dV 



where 



x(r) 



:£f (r)^(s)n(r) + £/(r)n(r) — Ara(r) 



-^r, e F (r) = f[3*Mr)} 2/3 , 
s F {r) 2m 



(3) 
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and we have used the universal form for the free energy 
per particle F/N in the unitary regime: 



F _ E-TS 
N ~ N 



3 3 
-e F ip{x) = ~e F [£,{x) 



xa(x)], (5) 



where for a homogeneous system — 5E/3e F N , 

<t(x) = S/N is the entropy per particle and x — T/e F 
(see inset in Fig. [1]). The overall chemical potential A 
and the temperature T are constant throughout the sys- 
tem. The density profile will depend on the shape of the 
trap as dictated by 6Sl/5n(r) — 0, which results in: 



sn 

Sn(r) 



5(F - XN) 
Sn(r) 



= H(x(r))+U(r)-X. (6) 



At a given T and A, equations (|4]) and ([6]) completely de- 
termine the density profile n(r) (and consequently both 
E(T,N) and S(T,N)) in a given trap for a given to- 
tal particle number. The only experimental input we 
have used is the particle number, the trapping potential 
and the scattering length at B = 1200 G, taken from 
Ref. 211 ] . The potential was assumed to be an 'isotropic' 
Gaussian, although it is not entirely clear to us to what 
extent this is accurate, especially in the axial direction. 
We have approximated the properties of the atomic cloud 
at B = 840 G with those at unitarity (B = 834 G), where 
we have MC data. For B = 840 G and for the parameters 
of the Duke experiment [2l| one obtains l/k F a = —0.06, 
using data of Ref. 30], if the Fermi momentum corre- 
sponds to the central density of the cloud at T = 0. 

Our results for the entropy of the cloud and the den- 
sity profiles for several temperatures, are shown in Figs. 
[2] and [3] In all the figures the temperature is measured 
in natural units of £_f(0), corresponding to the actual 
central density of the cloud at that specific tempera- 
ture. In Refs. 13, the temperature is expressed in 
units of the Fermi energy at T = in a harmonic trap: 
e h F ° = M7(3A0 1/3 . It is clear from Fig. [3]that the central 
density decreases with T and that the superfluid core 
disappears at T c = 0.23(2)£j?(0), which translates into 
T c = 0.27(3)4° to be compared to T c = 0.29(2)4° of 
Ref. [2l[ . There is a noticeable systematic difference be- 
tween theory and experiment at high energies, see Fig. [5] 
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FIG. 2: Entropy as a function of energy for the UFG in 
the Duke trap [21|: experiment (points with error bars) and 
present work (solid curve), where Eq = Ne 1 ^. Inset: loglog 
plot of E(T) as results from our calculations and as derived 
from experimental data The temperature is in natu- 

ral units, namely the Fermi energy at the center of the trap: 
£f(0). 
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FIG. 3: The radial (along shortest axis) density profiles of the 
Duke cloud at various temperatures, as determined theoreti- 
cally in the LDA using the MC results of Ref. The dot- 
ted blue line shows the superfluid part of the cloud, for which 
x(r) = T/ef{y) < 0.23. The solid red line shows the part of 
the system that is locally normal. Here a| = h/mui ma x- 



This discrepancy can be attributed to the fact that the 
experiment was performed slightly off resonance, on the 
BCS side, where l/k F a = —0.06. Even though theory 
will soon be extended to this region of l/k F a, a proper 
normalization of theory vs. experiment demands experi- 
mental results exactly at resonance. 

In Fig. [4] we show our results for the rms radius, in the 
form of the ratio of the mean square axial cloud size (z 2 ) 
at k F a = —0.75 to its value at unitarity, as a function of 
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FIG. 4: Ratio of the mean square cloud size (z 2 ) at akF = 
-0.75 (B = 1200 G) to its value near unitarity (B = 840 G), 
as a function of the energy: experiment [2lJ (points with error 
bars) and the present theory (solid blue line) . 



to Fig. [T]) . We might safely assume that the energy pre- 
diction based on MC data Q should not differ notice- 
ably from ours. However, a T c = 0.152(7)£j? will result 
in E c w 0.16i?o, in noticeable disagreement with exper- 
imental findings ID, 21]. Similarly, Ref. Q determines 
S(T C )/N = 0.2(2), which leads (according to our calcu- 
lations for the Duke trap) to S c — 1.5, noticeably less 
than the experimental value. One should notice that an 
analysis [3ll ] of the damping of sound modes [12] is also 
consistent with S c — 2.7. 

We are thankful to B. Clancy, L. Luo and J.E. Thomas, 
for useful discussions and for providing their experimen- 
tal data. We also thank Michael McNeil Forbes for ex- 
tensive discussions, in particular for suggesting the use 
of the function 1it{z). Support is acknowledged from the 
DOE under grants DE-FG02-97ER41014 and DE-FG02- 
00ER41132, and from the Polish Ministry of Science. 



the energy relative to the ground state. This dependence 
illustrates the relation between the energies of the cloud 
(related in turn to the spatial profiles via the virial theo- 
rem) at two different values of the magnetic field, but at 
the same value of the entropy. The quality of the agree- 
ment between our theoretical calculations and the exper- 
imental data demonstrates the soundness of the entire 
procedure to determine the entropy and the temperature 
scale for the UFG. In experiments, one can determine 
the value E c of the energy at the transition temperature 
without knowing the value of the temperature itself, sim- 
ply by noticing the appearance of a kink in E vs. S or 
vs. "empirical" T [131 ]. Specifically, it was determined 
in Refs. [H, that E c = E{T C ) - E(0) » 0.41(5)-E , 
to be compared with what we find theoretically for such 
a system E c = 0.32£^o- Similarly, it was determined in 
Ref. [IH that S c = S(T C )/N w 2.7(2), to be compared 
with our result S c = 2.15. We have identified T c with 
the disappearance of the superfluid core, which occurs 
according to our MC data at T c = 0.23(2)ej?, see Fig. 
|3J Our MC results are also consistent with a slightly 
higher T c » 0.25e_p, which would lead to E c = 0.36£'o 
and S c — 2.6 and thus to an almost "perfect" agreement 
between theory and experiment. As mentioned above 
however, an experiment exactly at unitarity is highly de- 
sirable, along with more precise MC data, in order to 
definitely settle the remaining discrepancies. The val- 
ues for the energy E c and entropy S c are particularly 
interesting because their determination does not require 
knowledge of T c and can be directly confronted with the 
MC calculations of Burovski et al. 0]. At low temper- 
atures, the lowest three values of E(T) of Ref. [§] agree 
with our own MC results The three highest tem- 
perature points, however, have in our opinion significant 
systematic errors (see Ref. [8] and comments in regards 
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